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Abstract: In this paper the operation of an amperometric biosensor producing a chemically 
amplified signal is modelled numerically. The chemical amplification is achieved by using 
synergistic substrates. The model is based on non-stationary reaction-diffusion equations. 
The model involves three layers (compartments): a layer of enzyme solution entrapped 
on the electrode surface, a dialysis membrane covering the enzyme layer and an outer 
diffusion layer which is modelled by the Nernst approach. The equation system is solved 
numerically by using the finite difference technique. The biosensor response and sensitivity 
are investigated by altering the model parameters influencing the enzyme kinetics as well as 
the mass transport by diffusion. The biosensor action was analyzed with a special emphasis 
to the effect of the chemical amplification. The simulation results qualitatively explain and 
confirm the experimentally observed effect of the synergistic substrates conversion on the 
biosensor response. 
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1. Introduction 

A biosensor is an electronic measuring device designed for measurement of the concentration of 
some specific substance (analyte) in a solution. The device specificity for a particular substance is 
achieved by some biological material, usually an enzyme [1-3]. An amperometric biosensor assesses 
the concentration of the analyte through the measurement of the current on a working electrode [4,5]. 
Biosensors are widely used in various applications that require fast quantitative analysis [6-10]. 

Oxidases are the kind of enzymes commonly used in biosensors. In the case of oxidases, molecules 
of an enzyme are reduced to their intermediate state and then oxidised by electron acceptors [5,11]. The 
reduction process as a rule is a specific reaction, while the oxidation of the reduced enzyme by electron 
acceptors is less specific. Many electron acceptors can be employed in this process [12]. 

The kinetic analysis of the reaction of a mixture of low and high reactive electron acceptors shows 
that the reduction of a low reactive electron acceptor may substantially increase if the rate of the cross 
reaction between these substances is high enough. This phenomenon is employed in biosensors for a 
high sensitivity determination of synergistic substrates [13-15]. 

The understanding of the kinetic peculiarities of the biosensors is of crucial importance for their 
design and optimization. To improve the productivity as well as the efficiency of the biosensor design, 
a model of the biosensor should be built [16,17]. Starting from the seventies various mathematical 
models have been widely used as important tools to study and optimize analytical characteristics of 
actual biosensors [18-22]. A comprehensive reviews on the modelling of the amperometric biosensors 
have been presented in [23,24]. Biosensors based on synergistic substrates were modelled at steady-state 
conditions a few years ago [25,26]. Very recently, a laccase-based biosensor utilizing simultaneous 
substrates conversion was numerically modelled at steady-state as well as transient conditions [27]. 

The goal of this investigation was to develop a computational model for an effective simulation of 
amperometric biosensors utilizing the synergistic substrates conversion as well as to investigate the 
influence of the physical and kinetic parameters on the biosensor response. The developed model is 
based on non- stationary reaction-diffusion equations [28,29]. The mathematical model comprises three 
compartments, a layer of enzyme solution entrapped on the electrode surface, a dialysis membrane 
covering the enzyme layer and an outer diffusion layer. By changing input parameters the output results 
were numerically analyzed at transition and steady-state conditions. The computational simulation 
was carried out using the finite difference technique [24,29,30]. The simulation results qualitatively 
explained and confirmed the experimentally observed effect of the synergistic substrates conversion on 
the biosensor response [13]. 

2. Mathematical Model 

We consider the following kinetic scheme used for the determination of synergistic substrates [13]: 

fci 

E ox + R >E Ted + P (la) 

Ered + S 1 -^^E ox + P 1 (lb) 
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E red + S 2 >E ox + P 2 (lc) 

Si + P 2 -^*Pi + S 2 (Id) 

Pi >S 1 + n 1 e~ (le) 

P 2 >S 2 + n 2 e~ (If) 



where P ox . and E re d stand for oxidized and reduced enzyme forms, respectively, R is the enzyme reducer, 
Si and S 2 are the substrates, P, Pi and P 2 are the products of the reactions, k\, k 2 , k 3 , fc 4 are the reaction 
rate constants, and ni, n 2 are the numbers of electrons involved in a charge transfer at the electrode 
surface in reactions (le) and (If), respectively. 

The reactions (la) and (lb) may be considered as the reduction of a low reactive electron acceptor Si 
by R and catalyzed by the enzyme. In the reaction (la), the enzyme is reduced by the reducer R, and an 
intermediate form of enzyme (E red ) is formed. In the reaction (lb), the electron acceptor Si is reduced 
by the E red , the product Pi is formed and enzyme molecules regain their primary oxidized form E ox . 

When the high reactive electron acceptor 5*2 is introduced, two more reactions, (lc) and (Id), start. 
The reaction (lc) represents the enzyme oxidation by the high reactive electron acceptor S 2 . The 
reaction (Id) is called a cross reaction. In the latter reaction, the substrate S\ oxidizes the product 
P 2 and regenerates the substrate S 2 . This reaction also produces the product Pi. 

Reactions (le) and (If) are electrochemical reactions that take place on the biosensor electrode. Both 
reactions are assumed fast and irreversible due to the high electrode potential. 

[Fe(CN) 6 ] 3 ~ ion may be employed as the low reactive electron acceptor Si [13,15]. In this 
case, [Fe(CN) 6 ] 4 ~ ion would be the reaction product P x . Cation radicals of l-(N,N-dimethylamino)- 
4-(4-morpholino)benzene, 5,10-dihydro-5,10-dimethylphenazine, phenoxazine-10-propionic acid or 
phenothiazine-10-propionic acid may be used as high reactive electron acceptors (5*2). The reduced 
forms of those substances would be considered as the corresponding products (P 2 ) [13]. 

During the physical experiment [13], the substrate 5*2 was introduced into the bulk solution at the 
moment when the steady-state response to the substrate Si was established. The thickness of the 
diffusion layer was controlled by the rotating magnetic bar and it was estimated using the rotating 
disk electrode. 

2.1. Biosensor Structure 

The biosensor is considered as an electrode with a relatively thin layer of enzyme solution trapped 
on the surface of the electrode by applying a dialysis membrane [13]. The biosensor model involves 
four regions: the enzyme layer where the enzymatic and chemical reactions as well as the mass 
transport by diffusion take place, a dialysis membrane, a diffusion limiting region where only chemical 
reactions as well as the mass transport by diffusion take place, and a convective region where the analyte 
concentration is maintained constant, see Figure 1 . Assuming a homogeneous distribution of the enzyme 
in the enzyme layer of the uniform thickness and symmetrical geometry of the dialysis membrane leads 
to the mathematical model of the biosensor action defined in a one-dimensional-in-space domain [23,24]. 
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Figure 1. Principal structure of a biosensor. 
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Let us define d\, d 2 and d 3 as the thicknesses of the enzyme, dialysis membrane and diffusion layers, 
respectively. We will also need values representing the distances between the electrode surface and the 
boundaries of the regions. Let a\, a 2 and a 3 be the distances between the electrode surface and one of 
those boundaries shown in Figure 1 . 

The diffusion layer (a 2 < x < a 3 ) may be treated as the Nemst diffusion layer [29]. According to the 
Nernst approach a layer of thickness d 3 remains unchanged with time. It was assumed that away from it 
the solution is uniform in concentration. 



2.2. Governing Equations 

The governing equations for a chemical reaction network can be formulated by the law of mass 
action [11,21]. Coupling reactions in the enzyme layer with the one-dimensional- in-space diffusion, 
described by the Fick's second law, leads to the following equations of the reaction-diffusion type 

(0 < x < ai, t > 0): 



9c red 

dt R1 dx 



-k x e ox r x + k 2 e red s n + k 3 e red s 21 (2a) 
heoxn - k 2 e red s u - k 3 e red s 21 (2b) 



(2c) 



~0f = D Sii-g^ 2 - - k 2 e red s u - kiS U p 21 (2d) 
ds 21 d 2 s 21 

~Qf = S2l ~dx2~ ~ k 3 e redS21 + Mllto (2e) 

-— - = D P^~Q^T + k 2^redS\l + hSup 21 (2f) 

dp2i n d 2 p 2l 

= Dp 2 x + k 3 e red s 21 - k A s u p 2 i (2g) 

where x stands for space, t is time, e ox (x,t) and e red (x,t) correspond to concentrations of oxidized 
(E ox ) and reduced (E red ) enzyme, respectively; s n (x,t) and s 21 (x,t) (p n (x,t) and p 2 i(x,t)) are the 
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concentrations of substrates (products) in the enzyme layer; ri(x, t) is the reducer concentration in the 
enzyme layer, and D R i, D Sl i, D S2 i, D Pll , Dp 2 i are the diffusion coefficients of the corresponding 
substances defined by the subscript. In the definition of the concentrations and diffusion coefficients, 
here and later in this paper the last numeric subscript label denotes the region of the model, particularly, 
1 stands for the enzyme layer. The molecules of both enzyme forms, E ox and E re d, are considered as 
immobilized, and therefore there are no diffusion terms in the corresponding equations. 

The product P does not act as a reactant in any reaction, so its concentration is not used in any further 
calculations. Therefore, equation system (2) contains no equation for the product P. 

No enzyme molecules appear in the dialysis membrane as well as in the diffusion layer. Only the 
reaction (Id) as well as the mass transport by diffusion of the reducer, substrates and products take 
place in these two regions. The governing equations for both these layers are represented as follows 
(cij_i < x < cii, t > 0, i = 2, 3): 

dri d 2 Ti 

m = Dm dx^ (3a) 

dsii n d2sii i nu\ 
— = D Sli -^- - k 4 s u p 2i (3b) 

ds 2i _ d 2 s 2 i . 

= D S2 i-^- + hsup2i (3c) 

dpu d 2 pu 

= D Pli -g-j- + k i s li p 2 i (3d) 

dp2i n d 2 p 2i 

= Dp 2i -g-j- - k±s u P2i (3e) 

where % = 2 corresponds to the dialysis membrane, and i — 3 corresponds to the diffusion layer. 

2.3. Initial Conditions 

Let x = 0 represents the electrode surface, while x = a\, x = a 2 and x = a 3 represent the boundaries 
between the adjacent regions as described in Section 2. 1 and shown in Figure 1 . The biosensor operation 
starts when the reducer and the substrates appear in the bulk solution. This is used in the initial conditions 

(t = 0), 

n(x,0) = s a (x, 0) =pii(x,0) = 0, 0<x<ai, £ = 1,2 (4a) 
e red (x, 0) = 0, e ox (x, 0) = e 0 , 0 < x < a x (4b) 

r 2 (x,0) = s i2 (x,0) = p i2 (x,0) = 0, oi < x < a 2 , £ = 1,2 (4c) 
Pi3(x, 0) = 0, a 2 < x < a 3 , £ = 1,2 (4d) 

r 3 (x, 0) = s i3 (x, 0) = 0, a 2 < x < a 3 , £ = 1, 2 (4e) 
^(^3,0) = r 0 , si 3 (a 3 ,0) = sio (4f) 

s 2 3K0) = { °' Ts ^l (4g) 

where e 0 stands for the total concentration of the enzyme in the enzyme layer (eo = e ox (x,t) + 
e ret i(x,t),Vx,t : x E (0, a a ) ,t > 0), r 0 is the reducer concentration in the bulk solution, s 10 and s 20 
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stand for substrate concentrations in the bulk solution, T$ 2 is the moment of the substrate S 2 insertion 
into the bulk solution. 

2.4. Matching Conditions 

On the boundary between two regions having different diffusivities, matching conditions have to be 
defined (t > 0, % = 1, 2, m = 1, 2), 



D 



dr r , 



Rm 



D 



dx 



Sim 



D 



Pivn 



dx 

dpim 



dx 



-D_R,m+l 
DSi,m+l 



dr m+1 



dx 



dx 

dpi, m +l 



dx 



j Pimiflm^i) Pi, m+l (tfm; ^) 



(5a) 
(5b) 
(5c) 



where m = 1 corresponds to the boundary between the enzyme layer and the dialysis membrane, 
whereas m = 2 corresponds to the boundary between the dialysis membrane and the diffusion layer. 

These conditions mean that fluxes of the reducer, substrates and products through one region are equal 
to the corresponding fluxes entering the surface of the neighboring region. Concentrations of the reducer, 
substrates and products in one region versus the neighboring region are assumed to be equal. 

2.5. Boundary Conditions 

In the bulk solution, concentrations of the reducer, substrates and products remain constant (t > 0), 



r 3 (a 3 ,t) = r 0 

s 13( a 3> t) — S W 

0, t < T S2 
S20, t> T S2 

#3(03,*) =0, i = l,2 



523(03, t) 



(6a) 
(6b) 

(6c) 

(6d) 



The reaction products Pi and P 2 take part in the electrochemical reactions (le) and (If), respectively, 
at the electrode surface (x — 0). Rates of those reactions are so high that the concentrations of P\ and 
P 2 at the electrode surface are permanently reduced to zero (t > 0) [13], 



Pu(0,t) = 0 
P2i(0,*) = 0 



(7a) 
(7b) 



Since the reaction (le) produces as much Si as it consumes Pi, the concentration flux of Pi on the 
electrode surface is equal to the concentration flux of Si but in opposite direction. The same is also 
assumed for the reaction (If) and substances S 2 and P 2 . These relations are expressed by the following 
boundary conditions (t > 0): 



D Pl i 



11 



dp 
dx 



-D 



2 = 0 



dsn 
dx 



(8a) 



x=0 
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Dp 2 i 



dp 2 i 
dx 



-D S2 i 



x=0 



dS2l 

dx 



(8b) 



x=0 



The reducer R is electrode-inactive substance, thus its concentration flux on the electrode surface is 
equal to zero (t > 0), 

dri 



D 



dx 



0 



(9) 



2 = 0 



2.6. Biosensor Response 



The measured current is usually assumed as the response of an amperometric biosensor in physical 
experiments. When modelling the biosensor action, due to the direct proportionality of the current to 
the area of the electrode surface, the current is often normalized with that area [23,31]. The biosensor 
current density j(t) at time t was expressed explicitly from Faraday's and Fick's laws, 



j 1 {t) = n 1 FD Pll 

h(t) = n 2 FD P2l — 
j(t)=ji(t)+j 2 (t) 



dx 
dp 2 i 



x=0 



x=0 



(10a) 

(10b) 
(10c) 



where ji(t) and j2(t) are the faradaic current densities generated by the electrochemical reactions (le) 
and (If), respectively, F is the Faraday constant, F = 96,486C/mol. 
We assume that the system approaches a steady state as t — > oo, 



lim j(t) 



t— >oo 



J St 

jut 

t— >oo 

j 2s t = Jim j 2 (t) 



lim ji(t) 



(11a) 
(lib) 
(He) 



where j st , j lst and j 2 st are the steady-state biosensor current densities. 

Since the current is generated due to two electrochemical reactions (le) and (If), it is important 
to investigate the influence of each of them to the overall biosensor response. The contribution of the 
electrochemical reaction (le) into the overall biosensor response was expressed as the ratio of the density 
of the steady-state current ji st to the density j st of the overall steady-state current, 



■h 



Jlst 
Jst 



(12) 



The sensitivity is also one of the most important characteristics of biosensors [1,2]. The biosensor 
sensitivity is usually expressed as the gradient of the biosensor current with respect to the concentration 
of the substrate in the bulk. Since the biosensor current as well as the substrate concentration varies even 
in orders of magnitude, a dimensionless expression of the sensitivity is preferable [31,32]. 

The biosensor based on synergistic reactions is designed to measure concentration of the substrate 
62 [13]. The dimensionless biosensor sensitivity B to the concentration of the substrate S2 was defined 
as follows: 



B(s 20 ) = x 



ds 



20 



Jst{S20) 
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where j s t(s 2 o) is the density of the steady-state biosensor current calculated at the concentration s 2 o of 
the substrate S 2 in the bulk. B(s 2 o) denotes the biosensor sensitivity at that concentration. 

The chemical signal amplification caused by a synergistic effect is another important characteristic 
of this type of biosensors [13]. A measure Amp of the synergistic effect was expressed as the ratio 
of the biosensor response to a substrate ^-containing analyte to the response to the corresponding S 2 - 
free analyte, 

Amp(s 20 ) = ^ (14) 

where Amp(s 20 ) is the dimensionless ratio of the signal amplification obtained by the insertion of the 
substrate S 2 of the concentration s 2 o into the buffer solution [13,22]. 

3. Numerical Simulation of Biosensor Action 

3.1. Simulating the Biosensor Operation 

Analytical solutions are not usually possible when solving non-linear partial differential 
equations [19,23,29]. Therefore, the biosensor action was simulated numerically. The simulation was 
carried out using the finite difference technique [29,30]. An explicit finite difference scheme was built 
on a uniform discrete grid with 200 points in space direction for each modelled layer corresponding to a 
certain time moment. The simulator has been programmed by authors in C programming language [33]. 

In the numerical simulation, the biosensor response time was assumed as the time when the change 
of the biosensor current remains very small during a relatively long term. A special dimensionless decay 
rate e was used, 

\dj(t) 



t 

t r = mm < t : 



m>o { ' j(t) 



<e}, j(t r )^j st (15) 



dt 

where t r is the biosensor response time. The decay rate value e = 1CT 3 was used in calculations. 

The following values of the model parameters were constant in all numerical simulations, unless 
stated otherwise [13]: 



fcl 


= 1.25 x 10 4 M" 1 




= 1.2 x 10 2 M~V 


h 


= 1.4 x 10 6 M -1 s 




1.4 x K^M^V 1 


e 0 


= 0.04mM, s w 


= 8mM, 


r 0 = 40 mM 


D R1 


= D Sil = D Pil = 


3.15 x 10" 


6 cm 2 /s, i = l,2 


D R2 


= D Sz2 = D Pi2 = 


4.2 x 1(T 7 


cm 2 /s, i — 1, 2 


D R3 


= D Sz3 = D Pi3 = 


6.3 x 1(T 6 


cm 2 /s, i — 1, 2 


di 


= 23.3 |^m, d 2 = 


= 18.6 |a.m 




m 


= n 2 = 1 







3.2. Experimental Validation 



The numerical solution of the model (2)-(10) was compared with the experimental data of the research 
conducted by one of the authors [13]. The results are depicted in Figure 2. The experimental data were 
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shifted by —30 s in time axis and by —1.9 (J.A in current axis [13]. The correction in time axis was applied 
in order to coincide the addition of reducer R (glucose in this particular case) into the bulk solution at the 
moment t — 0. The correction in current axis was applied in order to eliminate the background current 
influence to the overall biosensor response. 

Figure 2. Model validation with experimental data; (1) experimental data, (2)-(4) simulation 
results. Simulation was carried out at two reaction rate constants (fc 3 ) of the reaction (lc) 
(1.4 x 10 6 M" 1 s- 1 (2) and 0.35 x 10 6 M" 1 s" 1 (3, 4)) and at two thicknesses (d 3 ) of the 
diffusion layer (223 ^m (2, 3) and 116.5 ^m (4)); s 20 = 39 \iM. 

^ 20 i — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — 




400 500 



During physical experiment, the substrate S 2 was added into the bulk solution at the moment t = 130 s 
which corresponds to t = 100 s in the simulations. Thus T$ 2 = 100 s value was used in the numerical 
simulations. Having simulated density j of the biosensor current, the biosensor current was calculated 
by multiplying the density by the area A of the electrode surface, 

Ja = jA (16) 

In the physical experiment, the electrode with the surface area of A = 0.3 cm 2 was used [13]. 

The curve 2 in Figure 2 corresponds to the simulation carried out at the parameters equal to the 
parameters of the physical experiment [13]. The relative difference between the numerical solution 
and the experiment data at the moment T$ 2 (before insertion of S2) is about 14%. This indicates that 
mathematical model quite accurately reflects the physical experiment. However, after the insertion of 
S2, the accuracy deteriorates. The simulated steady-state current is about 2.5 times higher than the 
experimental one (Figure 2). This indicates that experiment is really more complex than as defined 
by the mathematical model. Presumably, the reaction of the reduced enzyme with oxygen, which 
was unforeseen in the mathematical modelling, is a main factor determining this difference [13]. The 
difference may be also caused by the limited stability of oxidized heterocyclic compounds and some 
external diffusion limitation of substrates [13]. 
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The curves 3 and 4 of Figure 2 demonstrate that simulation results may be fitted with the experimental 
ones by varying the value of the reaction rate constant k 3 and the diffusion layer thickness d 3 . However 
the measured values of the parameters will be used in this paper in further numerical simulations. 

Despite this inadequacy with the experiment, the model (2)-(10) seems to be suitable for investigating 
the kinetic peculiarities and optimizing the configuration of biosensors utilizing the synergistic effect 
of substrates. 

3.3. Simulating the Synergistic Effect 

The understanding the nature of the synergistic effect is of crucial importance for preparation of highly 
sensitive biosensors [13,15,26]. For a kinetic analysis of the biosensor operation, let us introduce V\, v 2 , 
v 3 , v A as the average rates of the reactions (la), (lb), (lc), (Id) taking place in the enzyme layer, 

vi = kie ox f (17a) 

v 2 = k 2 e Ted si (17b) 

^3 = k 3 e red s 2 (17c) 

v A = k 4 sip 2 (17d) 

where f, si, s 2 , pi, p 2 , e ox and e re d stand for the average concentrations of R, Si, S 2 , Pi, Pi, Eox and 
E red in the enzyme layer. 

Figures 3 and 4 show how the synergistic effect occurs. The biosensor parameters were kept the same 
as in Figure 2. 

Figure 3. Average concentrations f (1), s\ (2), s 2 (3), pi (4), p 2 (5), e ox (6) and e re d (7). 
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Figure 4. Average rates v\ (1), v 2 (2), v 3 (3), t> 4 (4) of the reactions (la)-(ld) taking place in 
the enzyme layer (curves 3 and 4 coincide). 

Q p 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 : 

= 1 

- 2 




0 0001 1 1 1 1 1 1 1 1 L ~' -1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

0 50 100 150 200 250 

t, s 

When S 2 is introduced into the bulk solution, a fast reaction (lc) starts, and a large amount of E ox and 
P 2 molecules are produced (see Figure 3). In the meantime the increased concentration of E ox initiates 
an increase in the rate of the reaction (la), molecules of P 2 start to react with Si in reaction (Id) (see 
Figure 4). A decrease of the reducer R concentration is the second indication of the increase in the 
reaction (la) rate as the reaction consumes more reactant molecules (see Figure 3). The reactions (la) 
and (lc) are very strongly interrelated because E ox is the product of reaction (lc), and at the same time 
it is the reactant of the reaction (la) and vice versa with E red molecule. Similar situation appears with 
the reactions (lc) and (Id) (molecules of P 2 and S 2 ). 

The rate of the reaction (lc) is controlled by the concentration of S 2 . Reactions (la) and (Id) are 
accelerated by the reaction (lc) and their rates tend to reach the rate of the reaction (lc). Hence the rates 
of reactions (la), (lc) and (Id) are approximately equal during almost all the time of biosensor operation 
as seen in Figure 4. The rate of the reaction (lb) remains almost unaffected. 

The reactions (lb) and (Id) produce the product Pi which is responsible for the biosensor current. 
Since the rate of the reaction (lb) remains almost unchanged and the rate of the reaction (Id) is greater 
than the rate of reaction (lb) by two orders of magnitude, the biosensor response is also increased by two 
orders of magnitude. In other words, the signal is amplified chemically. This is the desirable outcome of 
synergistic effect in this type of biosensors [13]. 

If the synergistic effect would not occur during the biosensor operation, that would mean that 
biosensor signal is not chemically amplified. A biosensor with no chemical amplification could be too 
weak for reliable measurement. That could mean the narrower range of biosensor application because 
biosensor would only work at relatively high concentrations of substrate. So, it is important to choose 
the correct biosensor parameters that support the occurrence of the synergistic effect. 



Sensors 2012, 12 



4908 



The law of mass conservation was applied for an additional validation of the numerical solution of 
the initial boundary value problem Equations (2)-(9). Since no mass transport of the enzyme takes place 
in the enzyme layer, the equality e ox (x, t) + e red (x, t) = e 0 holds for Va; G [0, ai] and Vt > 0. At the 
steady-state conditions when ds im /dt= dp im /dt = 0, the equality s im (x, t) + p im (x, t) = s i0 holds for 
Vx G [a m _i, a m ) and t — Y oo, i — 1, 2, m — 1, 2, 3. A fulfillment of these conditions can be observed in 
Figure 3 (please note that concentrations are provided in a logarithmic scale). 

3.4. Dynamics of Biosensor Current 

Figure 5 shows the dynamics of the biosensor current Ja simulated at two thicknesses (<i 3 ) of the 
diffusion layer and three concentrations («2o) of the substrate S2. For simplicity, here and below, the 
substrate S 2 is introduced in the bulk solution in the beginning of the experiments (Tg 2 = 0 s). 

Figure 5. The evolution of the biosensor current ja at different thicknesses (<i 3 ) of the 
diffusion layer (116.5 (xm (1, 2, 3) and 233 pm (4, 5, 6)) and concentrations (s 2 o) of the 
substrate S 2 (4 yiM (1, 4), 40 pM (2, 5) and 400 pM (3, 6). 




At low substrate concentration (s 2 o = 4 \xM) the steady-state current of the biosensor is not influenced 
by the thickness g? 3 of the diffusion layer. However at this concentration, a change in the thickness of 
the diffusion layer slightly changes the curvature of "current vs. time" curve and the biosensor with 
thinner diffusion layer reaches the steady state earlier. This is also correct in cases of higher substrate 
concentrations (S20 = 40 pM, S20 = 400 \iM). 

From the curves depicted in Figure 5 one can also observe another very important relationship: the 
steady-state current directly depends on the concentration s 2 o of the substrate S 2 in the bulk solution. 
At higher substrate concentrations (S20 = 40 pM, S20 = 400 pM) it is also observable that steady-state 
current depends inversely on the thickness d 3 of the diffusion layer. Intensity of the solution stirring 
plays important role when working with solutions of high concentration of the substrate. 
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4. Results and Discussion 



4.1. Biosensor Response vs. Substrates Concentrations 

The dependence of the steady-state biosensor current on the concentration S20 of the substrate S2 was 
investigated at different diffusion layer thicknesses d 3 as well as at different values of the reaction rate 
constant k 3 . The results of the numerical simulations are depicted in Figure 6. 

Figure 6. The dependence of the steady-state current density j st on the substrate concen- 
tration s 2 o at two thicknesses (d 3 ) of the diffusion layer (116.5 \xm (1, 2, 3) and 233 |xm 
(4, 5, 6)) and three reaction rate constants (k 3 ) of the reaction (lc) (7 x 10 5 M _1 s _1 (1, 4), 
1.4 x 10 6 M" 1 s- 1 (2, 5) and 2.8 x 10 6 M" 1 s" 1 (3, 6)). 




The biosensor may reliably measure the concentration of substrate S2 at the range of concentrations 
where the biosensor response is dependent on the concentration s 2 o but not dependent on the diffusion 
layer thickness d 3 . As one can observe from Figure 6 this range depends on the reaction rate constant 
k 3 . The lower boundary of this range varies from S20 ~ 0.02 |xM when k 3 is relatively high (curves 3, 6) 
to S20 ~ 0.1 \iM when k 3 is relatively low (curves 1,4). The upper boundary of this range varies from 
s 2 o ~ 16 (J.M when k 3 is relatively high (curves 3, 6) to s 2 o ~ 40 \±M when fc 3 is relatively low (curves 1, 
4). At the concentrations S20 > 250 |J.M the biosensor response is not influenced by the concentration of 
the substrate S 2 . 

To investigate the dependence of the biosensor sensitivity on the substrate concentration, the biosensor 
response was simulated at the same values of the diffusion layer thickness d 3 and the reaction rate 
constant k 3 as in Figure 6. The calculation results are depicted in Figure 7. 
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Figure 7. The dependence of the biosensor sensitivity B on the substrate concentration s 2 o- 
Values of all parameters are the same as in Figure 6. 




0.001 1 — — — — — 1 
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Figure 8. The dependence of the steady-state current density j st on the substrate 
concentration siq at two thicknesses (d 3 ) of the diffusion layer (116.5 urn (1, 2, 3) and 
233 \xm (4, 5, 6)) and three concentrations (s 2 o) of the substrate S 2 (1 M-M (1, 4), 10 uM (2, 
5) and 100 jiM (3, 6)). 




As one can observe from Figure 7, the biosensor shows the highest values of the sensitivity when the 
concentration of the substrate S 2 is moderate (1 |j.M < s 2 o < 10 M-M). At low substrate S 2 concentrations 
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(0.01 [iM < s 20 < 1 M-M) as well as at high concentrations of S 2 (10 \iM < s 20 < 1, 000 \xM) biosensor 
exhibits poor sensitivity. 

At low substrate concentrations, a higher k 3 yields a higher biosensor sensitivity, whereas the 
thickness d 3 of the diffusion layer does not influence the biosensor response at this concentration range. 
At high substrate concentrations, k 3 has the opposite effect on the biosensor sensitivity compared with a 
low substrate concentration range. A higher k 3 yields a lower biosensor sensitivity. At this concentration 
range, the thickness d 3 of the diffusion layer does influence the biosensor sensitivity. A thicker diffusion 
layer diminishes the value of biosensor sensitivity. 

Although this type biosensors are designed to measure the concentration s 2 o of the substrate S 2 , it is 
also important to investigate the biosensor response dependency on the concentration si 0 of the substrate 
Si in order to optimize the biosensors. The dependence of the steady-state biosensor current j st on the 
concentration s 10 was investigated at different values of diffusion layer thickness d 3 as well as at different 
concentrations (s 2 o) of the substrate S 2 ■ The simulation results are depicted in Figure 8. 

As we can observe from the curves in Figure 8, the biosensor response directly depends on the 
concentration si 0 . In cases of low and moderate substrate concentrations of S 2 , this dependency 
holds through the whole range of investigated s w concentrations. At high concentrations of £2 
(s 2 o = 100 (J.M), this dependency is valid only for low concentrations of Si (s w < 1 mM). At high 
concentrations of Si (s w > 1 mM) the biosensor response becomes practically insensitive to the change 
of s w . 

The thickness d 3 of the diffusion layer has no noticeable influence on the biosensor response at 
the majority of the biosensor parameters investigated. Only at high concentrations of Si and high 
concentrations of S 2 the effect of the diffusion layer thickness becomes noticeable. 

4.2. The Composition of Biosensor Response 

The biosensor response is determined by two electrochemical reactions (le) and (If). It is important 
to understand the role of each process at certain values of the biosensor parameters. 

The dependence of the biosensor response composition on the reaction rate constant k 4 was 
investigated at different diffusion layer thicknesses (d 3 ) as well as at different values of the substrate 
concentration s 2 o- Results of the numerical simulation are depicted in Figure 9. 

As one can see in Figure 9, the ratio J x of the density ji st of the steady-state current generated by the 
electrochemical reactions (le) to the density j st of the overall steady-state current directly depends on 
the reaction rate constant A; 4 . This can be explained by the reaction (Id). The faster the reaction (Id) the 
more P 2 consumes, and molecules of P 2 do not reach electrode surface. 

At low S20 concentrations, a value of £4 only slightly influences the composition of the biosensor 
response (curves 1 and 4). The biosensor response is almost entirely generated by the faradaic 
process (le) at such concentrations of S 2 . As one can see from the reaction scheme (1), a low 
concentration of S 2 directly cause a low concentration of P 2 . However at a high concentration s 2 o of 
S 2 , a value of the reaction rate constant A; 4 makes a major impact on the composition of the biosensor 
response. When fc 4 is small enough, the ratio Ji is less than 0.5, i.e., the main part of the biosensor 
response is generated by the faradaic process (If). 
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Figure 9. The dependence of the ratio Jjon the reaction rate k 4 at two thicknesses (d 3 ) of 
the diffusion layer (116.5 urn (1, 2, 3) and 233 \im (4, 5, 6)) and three concentrations (s 20 ) 
of the substrate S 2 (0.1 |oM (1, 4), 10 \±M (2, 5) and 1, 000 \iM (3, 6)). 




k 4 , M'V 1 

Figure 9 also shows that the thickness d 3 of the diffusion layer only slightly influences the composition 
of the biosensor response. 

4.3. The Role of Substrates Synergistic Conversion in the Biosensor Response 

The synergistic effect was analyzed with respect to the chemical amplification Amp. The dependence 
of the chemical amplification on the reaction rate constant k 3 was investigated at different diffusion layer 
thicknesses (e^) as well as at different values of the substrate concentration s 20 . Results of the numerical 
simulation are depicted in Figure 10. 

The direct influence of the reaction rate constant k 3 on the Amp is the main dependency that can be 
observed in Figure 10. At a higher rate constant (k 3 ) of the biochemical reaction (lc) more molecules of 
P 2 are produced. Because of this, the rate of the reaction (Id) increases, and eventually the amount of 
Pi also increases. Product Pi is the main generator of the current as was shown in the Section 4.2. 

As one can also see in Figure 10, a higher concentration s 20 of S 2 corresponds to a higher chemical 
amplification Amp. At low values of k 3 (k 3 < 10 3 Mr 1 s _1 ) this effect is rather weak, whereas at 
moderate values of k 3 (from 10 4 to 10 6 M _1 s _1 ) this effect is strongly expressed. However, at high 
values of the reaction rate constant k 3 (k 3 > 10 7 M _1 s _1 ) this effect noticeably diminishes again. A 
positive influence of s 2 o ° n me amplification can be explained by the increase of the rate of reaction (lc) 
that is caused by the increase of the concentration of one of the reactants. At low values of k 3 this is not 
so obvious because the increase in the reactant concentration does not counterweigh low values of k 3 . At 
high values of k 3 the system arrives at the point where the chemical amplification Amp is not dependent 
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on the reaction rate constant k 3 . The influence of the diffusion layer thickness d 3 appears noticeable only 
at high values of k 3 . 

Figure 10. The dependence of the chemical amplification Amp on the reaction rate constant 
k 3 at two thicknesses (d 3 ) of the diffusion layer (116.5 pm (1, 2, 3) and 233 pm (4, 5, 6)) and 
three concentrations (s 20 ) of S 2 (1 uM (1, 4), 10 [iM (2, 5) and 100 (3, 6)). 




Similar limitations of the rate of the synergistic process have been also observed experimentally using 
glucose oxidase biosensors utilizing the synergistic substrates conversion [13]. 

The amplification of the biosensor current was also investigated with respect to the substrate Si 
concentration. Numerical experiments were conducted at different diffusion layer thicknesses (d 3 ) as 
well as at different concentrations (s 2 o). The results of the numerical simulation are depicted in Figure 11. 

As one can see in Figure 11 the biosensor exhibits the highest values of the chemical signal 
amplification Amp at low concentrations of substrate Si and the lowest values of Amp at high values of 
Sio. So, the influence of sio on the biosensor response that was observed in Figure 8 is diminished by the 
chemical amplification, because the magnitude of the amplification buffers the impact of the substrate 
Si concentration change. When the concentration s w is low, the response of the biosensor is highly 
amplified by the synergistic effect and vice versa in the case of high sio- 

The thickness of the diffusion layer only slightly influences the chemical amplification throughout 
the range of the concentration s 10 that was used in this investigation. A slight influence of the diffusion 
layer thickness d 3 is observed only in the case of high concentrations of both substrates, Si and 5*2. 
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Figure 11. The dependence of the chemical amplification Amp on the substrate 
concentration s\ 0 at two thicknesses (d 3 ) of the diffusion layer (116.5 urn (1, 2, 3) and 
233 \xm (4, 5, 6)) and three concentrations (s 2 o) of the substrate S 2 (1 (J-M (1, 4), 10 uM 
(2, 5) and 100 uM (3, 6)). 




5 10 , raM 



5. Conclusions 

The developed mathematical model (2)-(10) is suitable for the modelling of an amperometric 
biosensor with synergistic substrates conversion. The model can be used as a tool for the biosensor 
optimization prior to the development stage of an actual biosensor. 

The investigation of the biosensor response showed that the range of concentrations, where the 
biosensor reliably operates, depends on the reaction rate constant k 3 . In cases when k 3 is lower the 
biosensor may be used to measure higher concentrations s 20 . In cases when k 3 is higher the biosensor 
may be used to measure lower concentrations s 20 (Figure 6). 

The biosensor has the highest sensitivity at moderate concentrations of the substrate S 2 (Figure 7). At 
low and high concentrations of S 2 , the biosensor sensitivity is quite low. 

The biosensor response is directly dependent on the concentration sio of the substrate Si. Low and 
moderate concentrations of Si are preferable, because at these concentrations the biosensor response is 
practically independent from the diffusion layer thickness (Figure 8). 

The biosensor response almost entirely depends on the electrochemical reaction (le) (Figure 9). The 
electrochemical reaction (If) may generate a significant current only in cases when the reaction rate 
constant of cross reaction (Id) is low and the concentration s 2 o of S 2 is moderate or high (Figure 9). 

The biosensor exhibits the highest values of the chemical amplification in cases when the reaction 
rate constant k 3 as well as the concentration s 20 are the high (Figure 10). The diffusion has a minor 
effect on the rate Amp of the chemical amplification. A thicker diffusion layer weakens the chemical 
amplification only at high values of the reaction rate constant k 3 and high concentrations of S 2 . The 
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amplification ratio Amp inversely depends on the concentration s w of the substrate Si (Figure 11). This 
effect is somewhat opposite to the si 0 effect on the biosensor response. 

To prove conclusions made the physical experiments are running using glucose dehydrogenase as 
well as laccase-based biosensors with different configurations. A more precise computational model 
implying the reaction of the reduced enzyme with oxygen in addition to the reactions (le) and (If) is 
now under development for a more accurate simulation of the biosensor action. 
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